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The average of the ratio of powers of the spectral determinants of the Dirac operator in the e-regime 
of QCD is shown to satisfy a Toda lattice equation. The quenched limit of this Toda lattice equation 
is obtained using the supersymmetric method. This super symmetric approach is then shown to be 
equivalent to taking the replica limit of the Toda lattice equation. Among other, the factorization of 
the microscopic spectral correlation functions of the QCD Dirac operator into fermionic and bosonic 
partition functions follows naturally from both approaches. While the replica approach relies on an 
analytic continuation in the number of flavors no such assumptions are made in the present approach 
where the numbers of flavors in the Toda lattice equation are strictly integer. 
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I. INTRODUCTION 



A quenched theory is one in which a determinant has been removed from weight in the partition function. Removing 
a determinant from the partition function does not necessarily lead to an unphysical theory. The determinant can 
act as a source term for Green's functions, and, once a derivative of the partition function has been taken, this source 
is removed. Two widely used analytical tools to remove determinants are the replica trick and the supersymmetric 
method. For the replica trick, one introduces n identical copies of the determinant, performs the appropriate deriva- 
tives, and finally takes the limit n — > 0. For the supersymmetric method, a ratio of two determinants with different 
values of the sources is introduced in the partition function. Derivatives are performed with respect to the sources 
and the determinants are removed by setting the sources equal. This paper deals with the equivalence of these two 
approaches. 

Originally, both the replica trick [1] and the supersymmetric method [2] where invented for the study of disordered 
systems. Applications of the replica trick soon proved to be delicate, see eg. [3], the trouble being in the choice of 
the dominant saddle points and the analytic continuation in the replica index [4]. The supersymmetric method [2] 
does not rely on such analytic continuation, and never faces this problem. At the perturbative level, one only finds a 
polynomial dependence on the number of replicas and both approaches yield identical results (see for example [4-6]). 
The replica trick has made its mark in diverse areas of theoretical physics. Some examples are the proof of the linked 
cluster theorem (see [7]) and the study of vacua in supersymmetric gauge theories [8]. 

In disordered systems, which also includes the Dirac operator in a gauge field distributed according to the Yang- 
Mills action, the standard approach both for the replica trick and the supersymmetric method is to rewrite the 
bosonic and/or fermionic integrals in terms of a non-linear cr-model. Such a cr-model describes the Goldstone modes 
associated with the spontaneous breaking of the global symmetries of the original model. Both its symmetry breaking 
term and its kinetic term are determined by the pattern of global symmetry breaking and Lorentz invariance. In the 
theory of disordered systems, these Goldstone modes are known as diffusons whereas in QCD they are known as the 
pseudoscalar mesons (for reviews see [2,9]). For QCD, the cr-model is known as a chiral Lagrangian and has been 
studied extensively in relation to the low-energy properties of the strong interactions [10]. 

In a parameter domain where the kinetic term is important, in most cases, only a perturbative treatment of the 
theory is feasible. However, in a domain where the kinetic term can be ignored, which is know as the ergodic domain or 
the e domain, the cr-model becomes zero dimensional and can be solved exactly in many cases. The zero dimensional 
a- models are equivalent to random matrix theories [2,11-14]. Using this connection, the problems with the replica 
trick where shown to be of a more serious nature [15]. This lead to the general consensus that the replica trick was 
only reliable for perturbative calculations and could generally not be used to obtain exact nonperturbative results. A 
first hint that nonperturbative results could be obtained from the replica trick came from [16,17] where the oscillatory 
factors in the asymptotic expansion of the two-point correlations function were reproduced. However, only in cases 
where the perturbative expansion terminates [18], exact analytical results could be obtained, and, as was pointed 
out in [18], a systematic and exact replica method remained illusive in these works. Recently, replica methods were 
introduced that went beyond these limitations [19-22]. Most significantly, exact analytical results were obtained in 
cases where other methods have failed. While the formulation of the replica trick introduced in [20] has been shown 
to be widely applicable, it still relies on an assumed analytical continuation in the replica index. The central purpose 
of this paper is to justify this analytical continuation. 

The new development which allowed for an exact analytic evaluation of the spectral correlation functions is the 
link between zero dimensional cr-models and an exactly solvable system know as the Toda lattice. This system is 
one dimensional and consists of massive particles interacting through an exponential potential [23]. The Hamilton 
equations of motion for the position of the masses as a function of time are known as the Toda lattice equations. The 
system is integrable, and its solution for appropriate boundary conditions is given by the partition function of a zero 
dimensional cr-model. 

Previously, we have shown [20,21] that the replica limit of the Toda lattice equation not only gives the exact analytic 
results but also explains the factorization of the universal microscopic correlation functions into a product of a bosonic 
and a fermionic partition function. Though this factorization suggests a direct relationship with the supersymmetric 
generating functional, it has remained a miracle why the final result factorizes within the supersymmetric method. 
In this paper we show this factorization based on explicit expressions for the ratio of different powers of spectral 
determinants [24]. The key observation is that this family of partition functions also satisfies Toda lattice equations. 
Quenching these Toda equations as in the supersymmetric method, that is setting bosonic and fermionic masses equal, 
automatically reveals the factorization of correlation functions. As stressed above, the replica limit of the Toda lattice 
equation relies on an analytic continuation in the number of flavors from integer values to real values. Although the 
replica limit of the Toda lattice equation, in cases that have been studied up to now, has given the correct analytic 
result, no definite arguments have been given that the usual problems of the replica limit do not manifest themselves 
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in this approach. As we will show in this paper the replica limit of the Toda lattice equation is equivalent to quenching 
the Toda lattice equation of a family of supersymmctric partition functions with a strictly integer number of bosonic 
and fermionic flavors. This result directly justifies the interpretation of the partition function with zero flavors in the 
usual Toda lattice equation as the supersymmetric partition function. 

Throughout this paper we consider the zero dimensional limit of the cr-model relevant for QCD. In the random 
matrix formulation this corresponds to the chiral unitary ensemble (chUE). We expect that similar results can be 
obtained for the unitary ensemble. 

This paper is organized as follows. First we describe the replica method and the supersymmetric method. Then, 
in section III, we discuss Toda lattice equations relevant for the one-point function and the two-point function. In 
section IV the quenching of the Toda lattice equations is discussed. This automatically yields the quenched correlation 
functions in the factorized form. The equivalence of the replica method and the supersymmetric approach is discussed 
in section V. Finally, in section VI, we summarize our results and comment on possible extensions of our work. In 
two appendices we derive the Toda lattice equations for the supersymmetric partition functions. 



II. DEFINITIONS OF REPLICA AND SUPERSYMMETRIC QUENCHING 

In this section we set our notation and introduce the supersymmetric method and the replica method. We illustrate 
the general definitions with generating functions for the microscopic limit of the one-point functions of the chUE. 

The essential ingredient in the partition functions are the fermionic and bosonic determinants which will be written 
explicitly whereas all other parts of the weight are included in the average denoted by (. . .). The fermionic replicated 
partition function is defined by 

Z Nf (x) = (det N t(D + x)), (1) 

where D is an operator, the mass a; is a number, and, as just mentioned, (. . .} denotes the average for the specific 
theory (note that d x (1) = 0). Analogously, the bosonic replicated partition function is given by 

Z- Nb {y) = (det- Nb {D + y)). (2) 

The replica trick for removing the determinants is to take the limit Nf — > or Nb — > 0. 

The supersymmetric partition function is defined by the average ratio of powers of the spectral determinants 

Z Nf ,- Nb (x\y) = (det N f(D + x)det- Nb (D + y)} (3) 

which can be removed by taking the limit y — > x. For example, 

lim Z 4 .- 2 (x\y) = Z 2 (x). (4) 

y— >x 

Obviously, it is not particularly interesting to introduce and remove determinants unless we differentiate with 
respect to the masses first. The average resolvent of the operator D is defined as 



In the fermionic replica limit it is given by 



G{x) = Jim Q -j^-d x \ogZ Nf (x), (6) 



in the bosonic replica limit we have that 



G{x) = lim ~krd x \ogZ- Nb {x), (7) 

N b ^0 — iV(, 



and the supersymmetric method is based on the identity 
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G{x) = lim0 x logZi,_i(z|y). (8) 

J/— >X 

While the average in eq. (1) appears to be well defined for any real value of Nf it is often not possible to evaluate 
the average unless Nf is an integer. For this reason, taking the limit Nf — > involves an analytic continuation in Nf. 

In this paper we focus on the topological trivial sector of the QCD partition function which is the average of a 
fermionic determinant 

ZN s {{m s }) = / [AAUY[ det(i 7 „D M - ro,)] e" 5 ™^ (9) 
J /=i 

over gauge field configurations with trivial topology weighted by the Yang-Mills action SVm- The quark masses 
{rrif} = mi, ■ ■ ■ , tunj may take different or identical values. In particular, we are interested in spectral correlation 
functions of the Dirac operator in the phase where chiral symmetry is spontaneously broken. The order 

parameter for the spontaneous breakdown of chiral symmetry is the chiral condensate denoted by (r/>^>). At low 
energies, the effective degrees of freedom are the Goldstone modes associated with this spontaneously broken symmetry, 
and the effective theory is a a- model [25]. The zero dimensional limit of this cr-model allows for an analytical evaluation 
of the correlation functions on microscopic scales. The microscopic scale is set by the inverse of Vfyip), where V is 
volume of the system. The domain in which the zero momentum modes of this effective theory decouple from the 
nonzero modes is known as the e-regime of QCD [10]. In the topologically trivial sector, v = 0, the fermionic partition 
functions are given by the unitary, and hence compact, integral [26-28] 

Z Nf ({x f })= [ dU e* Tc l M ' u+MU \ (10) 
Jueu(N f ) 

where M = diag(xi, • • • ,XN f ) = dia,g(('ip^p)Vmi 7 • • • , ($ip)Vm,N f ) is the rescaled mass matrix. The bosonic partition 
functions are given by a non-compact integral over positive definite matrices [13,29] 

Z-N b ({y b })= I dQe-i^Q+MQ-^ (n) 

JQeGl(N b )/U(N b ) 

where the dependence on the bosonic quark masses, the chiral condensate, and the volume is through the products 
ijb = mb{'tpip)V appearing in the rescaled quark mass matrix M = diag(yi, • • • ,yN b ). The supersymmetric partition 
function of QCD in the e-rcgime is given by an integral over a supergroup [12,13] 

ZNf,-N b (M\{y»})= I dQe* 8 ^^- 1 !, (12) 

JQeGl(Nf\N b ) 

where M — diag(a;i, ••• ,XN f ,yi, ■■■ ,yN b ) is the rescaled mass matrix. The boson-boson block of Gl(Nf\Nb) is 
Gl(Nb)/U(Nb) and the fermion-fermion block of Gl(Nf\Nb) is U(Nf). All three integrals are fixed uniquely by the 
symmetries. 

Much has been learned about these partition functions over the past few years. Explicit expressions in terms of 
determinants of Bessel functions are known in all cases [13,20,24,27,29,30,32-39]. For Z Nf and Z_ Nb the expressions 
have been derived directly from the group integrals (10) and (11) while the derivation of the general expression for 
ZN f ,-N b [24,34-36] goes through the random matrix representation. For all the partition functions Nf and Nf, must 
take integer values and the analytic continuation in Nf and Nb is not uniquely defined. 

As mentioned in the introduction, the latest developments include a link between such matrix integrals and an 
exactly solvable one dimensional system known as the Toda lattice [32,33,40]. The Hamiltonian equation of motion of 
the Toda lattice is known as the Toda lattice equation. The remarkable link is that the QCD partition functions also 
obey such Toda lattice equation. This equation connects partition functions with different numbers of flavors, much 
in the same way as recursion relations for orthogonal polynomials connect polynomials with different indices. It turns 
out that the Toda lattice equation for fermionic flavors is also valid for bosonic flavors. Both in the bosonic case and 
the fermionic case the Toda lattice equation terminates for zero flavors. Based on the correct analytical result for the 
replica limit of these hierarchies, it was conjectured in [20] that both hierarchies are connected by the supersymmetric 
partition function. This assertion is proved in this paper starting from the observation that the supersymmetric 
partition functions also obey a Toda lattice equation. 
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III. THE TODA LATTICE EQUATIONS FOR THE SUPERSYMMETRIC PARTITION FUNCTION 



The Toda lattice equations considered in the literature (see for instance [40-42]) connect partition functions with 
a different number of fermionic flavors or a different number of bosonic flavors. That is, the Toda lattice equations 
deal with semi-infinite hierarchies with either positive or negative index. In [20,21] it was found that these two 
semi-infinite hierarchies are connected. The connection was realized as the replica limit of the index in cither the 
positive or the negative semi-infinite hierarchy. Here, we consider Toda lattice equations for the partition functions 
with both fermions and bosons, i.e. super-symmetric partition functions. Referring to the graded symmetry of the 
supersymmetric partition functions we call these graded Toda lattice equations. In the next section we will derive the 
quenched resolvents from these graded Toda lattice equations. 



A. A graded Toda Lattice Equation Affecting One Flavor 

Starting from the expression for the partition function (12) with Nf fermionic flavors and N b bosonic flavors 
[34-36,20]* 

ZN S ,-N b {{xf}\{y b }) = — W - — — — T — -, (13) 

Uj>i=i( x j - x i) rij><=i(»f - VV 

where z t = x t for i = l,...,N f , z Nj+i = y, L for i = l,...,N b , Jj-i{zi) = ij-i(zj) for i = l,...,N f , and 
Jj-i(zN f +i) = {—^-) : '~ 1 Kj-i(yi) for i = l,...,N b , we prove in Appendix A that the partition function (12) sat- 
isfies the graded Toda lattice relation 

e xiv ft lir 1\ o 2 Z Nf + l,-N b {x l , {Xf}\{yb})ZN,-l,-N b ({Xf}fjii\{yb}) .... 

S Nf+Nb 5 Xi logZ N _ Nb ({x f }\{y b }) = 2nx, — / . (14) 

[ZN f -N b {{Xk}\{yb})\ 

Here, n is the number of fermionic flavors with mass Xi in ZN f ,-N b ({%f}\{yb})- The derivation 5 X is defined by 

S x = xd x (15) 

and 

f=l 6=1 

For degenerate masses the sum only extends over the different masses. 

This is the generalization of the Toda lattice equation for fermions obtained in [42] , to include both fermions and 
bosons. Below, in section IV A, wc quench this graded Toda lattice equation as in the supersymmetric method and 
show that the factorization of the quenched resolvent follows naturally. 



B. A graded Toda Lattice Equation Affecting Two Flavors 

The key ingredient in understanding the factorization property [21] of the two-point correlation function within the 
supersymmetric approach is a graded Toda lattice equation where two indices are raised and lowered simultaneously. 
In [21] it was shown that for m degenerate fermionic flavors with mass x\ and n degenerate fermionic flavors with 
mass X2, the partition function (10) satisfies the Toda lattice equation (m and n are positive integers) 



*Notation: The vertical bar in the argument of Z separates the masses {xf} = xi, . . . ,xn, of the fermionic quarks from the 
masses {y b } = j/i, . . . , j/jv 6 of the bosonic quarks. For degenerate masses we either repeat the mass in the argument of Z according 
to the degeneracy or put an additional subscript on Z indicating the degeneracy, eg. ^3,1,-2(^1, X2\y) = Zi-i(x\,x\,x\,x?\y, y)- 
For bosonic masses, the subscritpt of Z that denotes the number of flavors will contain and explicit minus sign. If all flavors 
are of the same kind the vertical bar will be omitted. 
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c c i ry { \ a 2 2 1 . n+ 1 1 ; #2 ) — 1 . n — 

l(Xl,X 2 ) 

^xi^ log Z min (a;i,a;2) = 4mnxf : ^- f — — i '-. (17) 

Using the explicit expressions for the partition function given in (13) we have shown that 

6 Xl 6 X2 l0gZ2-2(xi,X2\yi,y 2 ) = 4 Xl X 2 j- : rr^ (18) 

[Z 2 - 2 {xi,x 2 \yi,y 2 ))Y 

and a similar equation for other small numerical values of the number of flavors (see Appendix B). 

We conjecture that these results can be generalized to any integer number of fermions, Nf, and any integer number 
of bosons, Nb (further support for this conjecture will be given in Appendix B): 

x X , ry / 1 \\\ \\ a 2 2 Z N f +2,-N b ( x i, x 3A x f}\{yb})Z Nf -2,-N b ({ x f}f^iMyb}) nn , 

6 Xi 5 Xj log Z N _ Nb {{x f }\{y b }) =4nmx iXj =- , (19) 

[z Nf -N b {{ x f}\{yb})\ 

where m and n are the number of flavors with mass Xi and Xj in Z^ f ^N b ({xf}\{yi,}), respectively. A similar 
equation holds for differentiation with respect to two bosonic masses. In that case the fermionic flavors are passive. 
Furthermore, we have shown for a number of values of Nf and Nb (see Appendix B) that a corresponding equation 
is valid when we differentiate with respect to one fermionic and one bosonic mass. Without proof we assert that the 
general result in this case is given by 

x x ^ ^ (( ,,, -n , 2 2 z Nf+h-N b +i{x l7 {xf}\{y b }b^j)Z Nf - h - Nb - 1 ({xf}f^ l \y j ,{y b }) 

S Xi S Vj logZ Nf - Nb ({x f }\{y b }) = -Anmx iyj , (20) 

[ZNf,-N b ({Xf}\{yb})r 

where m and n are the number of flavors with mass Xi and yj in ZN f ^N b ({xf}\{yb}). These graded Toda lattice 
equations go beyond previous results [21] in two respects: 1) they have been extended to supersymmetric partition 
functions, and 2) the masses in these equations can be either degenerate or non degenerate. 



IV. SUPERSYMMETRIC QUENCHING OF THE GRADED TODA LATTICE EQUATION 

In the supersymmetric approach the partially quenched or fully quenched correlation functions are obtained by 
differentiating before taking the limit of equal fermionic and bosonic masses. In this section we show how the graded 
Toda lattice equations of the previous section automatically give the quenched one-point and two-point functions. 
The results are in exact agreement with expressions obtained previously. The important new result of this section is 
that the factorization of correlation functions into products of partition functions follows naturally without assuming 
an analytical continuation in the number of flavors. 



A. The One Point Function 



In this section we derive the quenched spectral one point function by quenching the graded Toda lattice equation 
(14) using the supersymmetric method. For this it is sufficient to have one bosonic and one fermionic flavor, and the 
graded Toda lattice equation then reads 



5x(&x + <y log^i ,-i(x\y) = 2x 



the limit y — > x of the l.h.s. of this equation is given by 



2 Z 2 , -i(x,x\y)Z- 1 (y) 



lim S X (5 X + 6 y )logZi- 1 (x\y) = lim 6 X 



x( Tr- 



1 



D + x 



y(Tr 



D + y 



■r,y 



Here, we have used the notation 



(• • -)x,y 



8 x xG(x). 

dct{D + x) 
' det(L> + y) 



(21) 



(22) 



(23) 



G 



Taking also the limit y — > x of the r.h.s. of (21) we find 

5 x xG{x) = 2x 2 Z 1 {x)Z- l {x). (24) 

This result is in complete agreement with known results for the resolvent and the partition functions. For completeness 
we mention that they are given by 

G{x)=x(K {x)h{x) + K_ 1 {x)h{x)), Z 1 {x)=I {x), Z_ x (x) = K (x). (25) 

Clearly, this is a rather cumbersome way of deriving the resolvent from Z\-\. Let us therefore again stress that the 
important new result here is that the factorization in eq. (24) is proved on the basis of the Toda lattice hierarchy for 
the supersymmetric partition functions. We also stress that this approach does not rely on an analytic continuation 
in the number of flavors. 

Following the same strategy as above one can also obtain the partially quenched resolvents. As we seek to stress a 
matter of principle rather than reproducing existing results we instead turn to the two-point function. 



B. The Two-Point Function 



The two-point disconnected susceptibility contains two derivatives just as in the Toda lattice equation. This property 
simplifies the derivation of two point-functions from the Toda lattice equation for ratios of spectral determinants. 
The fully quenched disconnected susceptibility in the chiral unitary ensemble with identical arguments is defined as 

X(x,x) = - lim d x d v logZ lt _ 1 (x\y). (26) 

y^x 

It follows from (20) for Nf = N b = 1 at equal values of the arguments that 

X{x,x) = - lim. d x d y logZi.-i(x\y) — Ax 2 Z 2 (x, x)Z- 2 (x, x). (27) 

y^x 

In order to get the quenched disconnected susceptibility for different values of the arguments, we start from the 
partition function describing two fermions and two bosons all with different masses. We now have a choice: we can 
differentiate with respect to two fermionic masses, with respect to two bosonic masses, or with respect to one fermionic 
mass and one bosonic mass, before pairing up the fermion and boson masses. If we choose the first option and use 
(18) we find 

X(xi,yi) = ~ lim d xl d X2 \ogZ 2 ^ 2 {xi,x 2 \y 1 ,y 2 ) (28) 

2/1— >a:i,J/2— >x 2 

yi^x 1 ,y 2 ^x 2 [Z 2 ,- 2 (xi,x 2 \yi,y 2 )] 2 
= A.X\X 2 Z 2 {x\,x 2 )Z_ 2 {x\,x 2 ). 

The derivation when differentiating with respect to the two bosonic masses is analogous and leads to exactly the same 
quenched correlation function. If we choose to differentiate with respect to one fermionic mass and one bosonic mass 
and use (20) for N f = N b = 2, 

o a 1 ry / | \ „ Z 3i _ 1 (x 1 ,x 2 ,x 2 \y 2 )Z 1 ^ 3 (x 1 \y 1 ,y 1 ,y 2 ) 

d X2 dy 1 logZ2-2{x 1 ,x 2 \yi,y2) = -4x 2 yi r- : \y2 ' ( 29 ) 

[Z 2 - 2 (xi,x 2 \yi,y 2 )\ A 

(which we have also verified explicitly) we obtain 

dx 2 

2 

,. , Z 3 - 1 (x 1 ,x 2 ,x 2 \y 2 )Z 1 .- 3 (x 1 \y 1 ,y 1 ,y 2 ) 

hm 4:x 2yi — — — 

yi^x uy2 ^x 2 [Z 2 - 2 (xi,x 2 \yi,y 2 )\ z 

= kx\x 2 Z 2 {x\,x 2 )Z_ 2 {x\,x 2 ). 

This result coincides with (28) and is in agreement with the result obtained from a supersymmetric calculation [43] 
and with the result obtained from the replica limit of the Toda lattice equation of the fermionic hierarchy [21]. Most 
importantly, it is clear that the factorization property of the quenched two-point function in the supersymmetric 
approach is a direct consequence of the structure in the graded Toda lattice equation. No assumptions on the analytic 
continuation in the flavor index have to be made in this case. 



X(xi,x 2 ) = - lim dx 2 d Vl logZ 2 - 2 (x 1 ,x 2 \y 1 ,y 2 ) (30) 

3/1— >Xl,3/2— >X2 
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V. EQUIVALENCE OF REPLICA AND SUPERSYMMETRIC QUENCHING 



In this section we show that the replica limit of the Toda lattice equation necessarily gives the same answer as the 
supersymmetrically quenched Toda lattice equation. This directly justifies the central assumption of [20,21], namely 
that the bosonic and fermionic semi-infinite Toda lattice hierarchies are connected by the supersymmetric partition 
function. 



A. Equivalence for the one Point Function 

The QCD partition function (10) with Nf mass degenerate fermionic quarks satisfies the Toda lattice equation [42] 

S x \ogZ Nf (x) = 2N f x ^ • (31) 

(Note that this is obtained from (14), satisfied by the supersymmetric partition function, setting n = Nf and N b = 0.) 
In the replica approach we would like to take the Nf — * limit of this equation. Using the representation (1) of the 
partition function, which depends continuously on Nf, we get [20] 

lim \ogZ Nf (x) = S x xG(x), (32) 

I\ f — >v) 1\ f 

lim 2x 2 ^p^zM = 2x 2 Zl{x)z _ l{x y (33) 
N f ^o [Z Nf {x)] 2 

However, the Toda lattice equation (31) was established starting from (10), which is only well defined for positive 
integer values of Nf, and the boundary condition Zjf f=0 (x) = 1. The Toda lattice equation (31) therefore is not 
necessarily valid for (1) when Nf is non-integer. However, as we will now show, equating (32) to (33) is identical to 
quenching the Toda lattice equation using the supersymmetric method, since the supersymmetric partition function 
satisfies the Toda lattice equation. To this end we write (31) as 

1 c 2n ry / \ 2 Z Nf +i{x)Z N x {x) Z N +2 ,-i(x\y)Z N ^(xly) 

— 8 Z X log Z N (x) = 2x A '—r- — = lim 2x z '—r- — '— ,2 . (34) 

Nf [Zn s (x)] 2 y^x [Z N!+1 ^{x\y)\ 2 

These equalities are valid for positive integer values of Nf. However, the r.h.s. of the second equality is also well 
defined for Nf = 0. If we take this as the definition of the Nf — > limit of the middle term of this equation we 
get 2x 2 Zi(x)Z^i(x) exactly as in (33), but now without ambiguities in the interpretation of Z_\(x). Moreover, in 
appendix A it is shown that the supersymmetric partition function entering the r.h.s satisfies the same Toda lattice 
equation 

X fX , X M 7 f \ \ OfAT , -n 2 Z N,+2,-l( X \y) Z N f ,-l(x\y) . . 

S X {S X +S y ) log Z Nf+ly -!(x\y) = 2(N f + l)x — n ^ . (35) 

[ZN f +i-i(x\y)\ 

This equation is also well-defined for Nf = and in this case the limit y — > x is (cf. (22)) 

lim 5 X (S X + 6 y ) log Zi i(x\y) = S x xG(x). (36) 

y—>x 

This is precisely what we obtained in (32) using the representation (1) of the fermionic partition function. This shows 
that the replica limit of the Toda lattice equation is necessarily identical to the quenched limit of the supersymmetric 
Toda lattice equation. 



B. Equivalence for the Two-Point Function 

We now consider the replica limit of the Toda lattice equation affecting two indices simultaneously. In [21] it was 
shown that for integer n > 

a a l [7 / M a 2 Z n+ i n+ i(xi,X2)Z n _i n _i(xi,X2) , . 

d xl d X2 \og[Z n , n (x 1 ,x 2 )\ = 4n z x 1 x 2 : — — — . (37) 

[Z n ,n{Xl,X2)\ Z 
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Making the analytic continuation in n it was assumed [21] that the limit n — > connects this semi-infinite hierarchy 
to the semi-infinite hierarchy for Z_ n> _ n , viz. 

lim \d Xl d X2 \og[Z nn {x l ,X2)\ = 4 : x 1 x 2 Z 1 . 1 (x 1 , x 2 )Z_ 1 -i(x 1 , x 2 ). (38) 

n^O n 

While the exact agreement of the resulting quenched two-point correlation function with the known result [43] vindi- 
cates this replica limit, the underlying assumption again could not be proved. 

We now prove the assumption for the two-point function using similar arguments as for the one-point function. Again 
we establish the troublesome connection between positive and negative integer values of n by quenching determinants 
as in the supersymmetric method. The ratio of partition functions in the r.h.s. of (37) can be written as 

Z n+lin+ i(xi,x 2 )Z n - ltn - 1 (x 1 ,x 2 ) _ Z n+2 ^ +2 - 1 _i(x 1 ,x 2 \y ll y 2 )Z n ^ n ^ 1 _i(x 1 ,x 2 \yi 1 y 2 ) 
[Z n ,n(Xl,X 2 )] 2 xi.l/2— x 2 [Z n+ i tn+ i-i-i(x 1 ,x 2 \yi,y 2 )} 2 

As discussed in section IV B the ratio on the r.h.s. satisfies a Toda recursion relation even before taking the limit of 
degenerate masses 

4(n+l) x x x 2 ; '-jj= : = d Xl d X2 log [Z n+1 , n+1 _i _ 1 (xi,x 2 3/1,3/2) (40) 

[Z n+ i )n+ i _i _i(xi, x 2 \yi, y 2 )\ 2 

In this relation we can safely replace n by zero so that 

Y\m \d Xl d X2 \og[Z ntn {x\,x 2 )]= lim d Xl d X2 log[Z"i,i _i ,-i(xi,x 2 \yi,y 2 )\. (41) 

n^O n J/i— >xi,j/2— >X2 

This shows that the replica limit of the Toda lattice equation necessarily gives the same answer as the supersymmet- 
rically quenched Toda lattice equation, thus justifying the analytic continuation in n. 



VI. CONCLUSION 



The QCD partition function in the e regime is equivalent to the microscopic limit of a chiral random matrix theory 
with the same global symmetries. This equivalence holds for any number of bosonic and fermionic quarks, and 
each partition function in this family can be expressed as a determinant of Bessel functions. In this paper we have 
shown that this family of partition functions satisfies a Toda lattice equation which we have called the graded Toda 
lattice equation. Based on this result we have shown that the semi-infinite hierarchy of bosonic or fermionic partition 
functions can be continued to the limit of zero flavors which is given by the supersymmetric partition function. This 
was the central assumption when taking the replica limit of the Toda lattice equation. By proving this assumption 
we have thus shown the equivalence of the replica limit of the Toda lattice equation and the supersymmetric method. 
Quenching the partition functions in the graded Toda lattice equation shows that the resolvent of the supersymmetric 
partition function factorizes into the product of a bosonic and a fermionic partition function. 

In its simplest form, the graded Toda lattice equation relates derivatives of the supersymmetric partition function 
to partition functions with one more flavor or with one flavor less. This case is relevant for the analysis of the one- 
point function and has been proved in all its generality. For the analysis of the two-point function we need a graded 
Toda lattice equation in which derivatives of the partition function are related to partition functions with two more 
flavors or two less flavors. In addition to the cases that were required for our analysis, this Toda lattice equation was 
proved for several other cases with a fixed number of flavors. Based on these results and known results for exclusively 
fermionic flavors, we have conjectured that this equation can be generalized to an arbitrary number of bosonic and 
fermionic flavors. 

It is not yet clear whether variants of the Toda lattice equation are suitable for the investigation of the replica 
limit of real symmetric ((3d = 1) and quaternion real ((3d = 4) random matrix ensembles. However, our work can be 
extended in several other directions. First, our results can be extended in a straightforward way to partition functions 
with a nonzero topological charge. The reason that we considered only zero topological charge in this paper is that it 
somewhat simplifies our notation. Second, in the formulation of Kanzieper, the replica limit is taken of the Painlcvc 
equation rather than the Toda lattice equations. It would be very interesting to understand if the supersymmetric 
partition functions can also be incorporated in this framework. Third, the spectral density of the quenched Dirac 
operator at non-zero chemical potential in the limit of weak non-hermiticity was obtained from the replica limit of a 
Toda lattice equation. We expect that the supersymmetric generating functional at non-zero chemical potential will 



9 



satisfy a similar Toda lattice equation so that the replica limit can also be justified in this case. Fourth, the two-point 
function of the unitary ensemble can also be obtained from the replica limit of a Toda lattice equation. Also in this 
case we believe that a generalization of the Toda lattice equation to supersymmetric partition functions is possible. 
Finally, we have no doubt that a similar arguments can be applied to QCD in 3 dimensions. 
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DE-FG-88ER40388. 



Appendix A. THE GRADED TODA LATTICE EQUATION WHERE ONE INDEX IS AFFECTED 

In this appendix we prove the graded Toda lattice equation (14) relevant for the partially quenched one point 
functions, that is the relation where one flavor is added and subtracted. We start by slightly rewriting the general 
expression for the supersymmetric partition function in (13) [34-36,20] 



v ,r xv \a dct[<^ 1 J(z i )] itj=h .. tNf+Nb 

IL>i=i(ay - x i) Y\j>i=M - Vi) 



(42) 



Here, Zi = Xi for i = l,...,Nf, and ZN f +i = y% for i = l,...,N b . Furthermore, the symbol J{z{) = Io(%i) for 
i = 1, . . . , Nf and J(zN f +i) = -Ko(&/i) for i = 1, . . . , N b . In order to keep track of the indices let us consider the last 
two fermionic flavors as mass degenerate, that is xn } — Xn s +i- In this case the partition function is equal to 



Z Nf+1 _ Nb (x Nf ,{xf}f= lt .. tNf \{y b }) = 



det[A Nf+1 _ Nb (x Nf ,{x f }f =lj .. jNf \{y b })} 



where we have introduced the (Nf + 1 + AT&) x (Nf + 1 + Nb) matrix 

/ Io{xi) S Xl I (xi) 



A Nf+1 _ Nb (x Nf ,{x f } f=l! .. iNf \{y b }) = 



Io(xN f ) 

d XNf Io(xN f ) 
K (Vi) 

V K o(yN b ) 



8yiKo(yi) 



c +JV %(*l) \ 

Sx^ Nb Io(xN f ) 

o XNf 6 x ' I (x Nf ) 



SyJ b +NbK o(yN b ) J 



(43) 



(44) 



The proof of the graded Toda lattice equation follows from the Sylvester identity [44] valid for the determinant of any 
matrix 



^i>j^p>q ^i>q^P;j — det (A^C^j-pq , 
where Cjj is the cofactor of matrix element ij 

ddet(A) 



dAi 



and Cij-p yq is the double cofactor of matrix elements ij and pq 

d 2 det (A) 



(45) 



(46) 



(47) 



Setting i = N f + 1, j = N f + 1 + N b , p = Nf, and q = N f + N b we have 

CN f + l,N f + l+N b CN f ,Nf+N b — CN f + l,Nf+N b CN f ,Nf + l+N b = det(Ajv / +l,-AT ! , )CN f + l,N f + l+N b ;Nf,N f +N b ■ (48) 
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From the structure (44) of the matrix i4jv / +i,jv b it follows that (SN f +N b is defined in (16)) 

CN f + l,N f + l+N b = detAN fi _N b , 

CN f .N f +N b = S Xn 8N f +N b dct AN f -N b , 

X Nj 

CN f +i,N f +N b = $N f +N b det A N; _ Nb , 

CN f ,N f + l+N b = S x det AN f ,-N b , 

XNt ' 

CNf + l,Nf + l+N b ;Nf,N f +N b = det AN f -l,-N b - 

For completeness we write out 

/ I {xi) S Xl I (xi) ■■■ 5 Xl f ~ 1+Nb I (x 1 ) \ 



A Nf ,-N b {{xf}\{yb}) = 



h {XN f ) 

K {yi) S Vl K ( yi ) 



x N f -l + N b T , , 



rN f -l + N b 



K (yi) 



S^r 1+Nb Ko(y Nb ) J 



\ K {y Nb ) 
Inserting this in (48) we get 

[det A Nf _ Nb ] 2 5 XNf S Nf+Nb log[detA Nf - Nb ] = x Nf Act A Nf+ i_ Nb det A Nf _i_ Nb 



(49) 



(50) 



(51) 



Finally, using the relation (43) between the partition functions and the matrices A we obtain the desired graded Toda 
lattice equation 



SN f +N b S Xt \ogZ Nf - Nb ({x f }\{y b }) = 2x 



2 Z Nf+ i_ Nb (x i ,{xf}\{y b })Z Nf _ 1 _ Nb ({xf} fiH \{y b }) 



[Z Nf ,-N b ({x k }\{y b }W 



(52) 



Note that the Vandermonde determinants, A({zf}), in the squared masses do not contribute to the left hand side of 
the Toda lattice equation since SN f+ N b log[A({^?})] is a number. 



In the case where n flavors have mass xjy f the proof is similar. The partition function is now given by 
, , , ., det[A n , N h - Nb (x Nf Axf}f^N f \{yb})} 

Zn,N f -l,-N b {XN f ,{ x f}f?Kf\{ybl) ~ 



i\%=m - x\) n?Ji~\x% f - x*)"-n2x Nt )«»-i)/2 nS=i(»i - yf) 

where we have introduced the (Nf + n — 1 + N b ) x (Nf + n — 1 + N b ) matrix 

/ Jo(xi) S Xl I ( Xl ) ■■■ S^ +n - 2+Nb I ( Xl ) \ 



, (53) 



A n ,N f -l,-N b (XN f ,{Xf}f^N f \{yb}) 



Io(x Nf ) 
jid XNf I (x Nf ) 

j^Ty^xN^oixNf) 

Ko(yi) 8 yi K (y 1 ) 



5^ n - 2+N *Io(x Nf ) 
±d XNj 5^ n - 2+Nh Io(x Nf ) 



1 an-ls:N f +n-2+N b j , v 
J^=TyO XNf XNf h( x Nj) 

6 N f+n -2 + N bKo{yi) 



x Nj+n-2+N b „ , s 



(54) 



V K (y Nb ) 

Using again the Sylvester identity, (45), this time with i — Nf + n — 1, j = Nf + n — 1 + N b , p — Nf + n — 2, and 
q = Nf + n — 2 + N b we find the graded Toda lattice equation 

(>N f +N b S Xi logZn^-i.-jVfcfai, {xf}f^i\{y b }) (55) 

2 Zn+l,N f -l,-N b { x i>{ x f}f¥ : i\{yb}) Z n-l,Nf-l,-N b { x u{ x f}f¥ : i\{yb}) 



2nxf 



[Zn,N f -l,-N b (Xi,{Xk}fjti\{yb}W 
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Appendix B. THE GRADED TODA LATTICE EQUATION WHERE TWO INDICES ARE AFFECTED 



In this appendix we consider the graded Toda lattice equations (19) and (20) relevant for the two-point functions. 
Although we have been able to check these equations in a number of cases including those relevant for quenching the 
Toda lattice equation we do not have a general proof. 

The simplest possible example is 

d x d v \o g z lt _My) = -^v ^mI)] 2 ' (56) 

This result follows by rewriting the l.h.s. of this equation as 

[z^Mv)?o x d y io gZl ,_ M y) - det ( fX-% IfX-Mv) ) ' (57) 



and using the identity 

Zi,-Mv) d x Z h _My) \ = lf Z i( x ) dxZi(x) \ [0 1 \f Z- X {y) d y Z^(y) 

d y Z lj _ 1 (x\y) d v d x Z x .^{x\v) j 2 \ d x Z^{x) d 2 x Z x (x) J { 1 J \ d y Z^(y) d 2 y Z^{y) 

which valid if the partition function is of the form 



(58) 



Zi,-i{x\y) = X -[Z x {x)5 y Z^{y) + b x Z x {x)Z_ x {y)\. (59) 

This structure of the supersymmetric partition function easily follows by integrating the bosonic variables before the 
Grassmann variables, 



Z\,-MV) = J dadf3I (x(l + a[3/2)K (y(l -a/3/2) 



-[I (x)6 y K (y) + 8 x I (x)K (y)]. (60) 



The general case: Our conjecture is that the generalization of the above graded Toda lattice equation is given by 

X r , ry (( Mf ^ A 2 2 Z Nf + 1 ,- N b + 1 fa, {x f } | {y b } h+ j ) Z N _ 1 , _ N>> _ 1 ( {x f } | J/j , {y b } ) 

S Xt 5 Vi \ogZ N f ,-N b ({x f } {y b }) = -4nmx tyj — '- — / , (61) 

where n is the degeneracy of Xi and m is the degeneracy of yj in Z N ft -N b ({xf}\{yb}), respectively. We have explicitly 
checked this equation for the cases listed in the table below: 



Eq (61) 


N f 


N b 


n 


m 


case 1 


1 


1 


1 


1 


case 2 


2 


2 


1 


1 


case 3 


2 


2 


2 


1 


case 4 


2 


2 


2 


2 


case 5 


3 


1 


3 


1 


case 6 


1 


3 


1 


3 



A similar equation holds when both derivatives are with respect to fermionic masses. It is conjectured to be given 

by 



x x , ry ff n, n , 2 2 Z N s +2,-N b {x l iX J1 {x f }\{y b })Z N f -2,-N b {{xf}f^i,j\{y b }) 

5 Xi 5 Xj logZ N _ Nb ({x f }\{y b }) = \nmx,x J — =- , (62) 

[ZN s -N b {{Xf\\{y b \)\ 



where n is the degeneracy of Xi and m is the degeneracy of Xj in Z^ f ,-N b ({xf}\{y b }), in this order. We have verified 
this identity for the following cases: 
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Eq (62) 


N f 




n 


m 


case 1 


any 





any 


any 


case 2 


o 
Z 


o 

z 


1 


1 


case 3 


3 


1 


2 


1 


case 4 


4 


1 


2 


2 


case 5 


4 


1 


3 


1 
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